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QCD sum rules analysis of the B s — > D s j(2460)Zz/ decay 
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Using three point QCD sum rules method, the form factors relevant to the semileptonic B s — > D a j(2460)£v 
decay are calculated. The q 2 dependencies of these form factors are evaluated. The dependence of the asymmetry 
parameter a, characterizing the polarization of D s j meson, on q 2 is studied. This study gives useful information 
about the structure of the D s j meson. Finally the branching ratio of this decay is also estimated and is shown 
that it can be easily detected at LHC. Note that this talk is based on the original work in 1]. 

I. INTRODUCTION 

Recently, very exciting experimental results have been obtained in charmed hadron spectroscopy. The obser- 
vation of two narrow resonances with charm and strangeness, D s j(2317) in the D s ir° invariant mass distribution 



mm. 



and -D s j(2460) in the D*w and D s ^ mass distribution 



sions about the nature of these states and their quark content 



10 



3,0, y,Q,y,3, 



has raised discus- 



Ill. Analysis of the D So (2317) -> D*j, 



Z? s ,/(2460) — » Z)*7 and I?5 J (2460) — > D So (2317)7 indicates that the quark content of these mesons are probably 



cs 



12]. In 



12| it is also shown that finite quark mass effects for the c-quark give non-negligible corrections. 
When LHC begins operation, an abundant number of B a mesons will be produced creating a real possibility 
for studying the properties of B s meson and its various decay channels. One of the possible decay channels of 
B s meson is its semileptonic B s — ► D s j(2460)/z/ decay. Analysis of this decay might yield useful information 
for understanding the structure of the D s j (2460) meson. 

It is well known that the semileptonic decays of heavy flavored mesons are very promising tools for the 
determination of the elements of the CKM matrix, leptonic decay constants as well as the origin of the CP 
violation. In semileptonic decays the long distance dynamics are parameterized by transition form factors, 
calculation of which is a central problem for these decays. 

Obviously, for the calculation of the transition form factors, nonperturbative approaches are needed. Among 



the nonperturbative approaches, the QCD sum rules method 



1 31 ] received special attention, because this method 



is based on the fundamental QCD Lagrangian. This method has been successfully applied to a wide variety of 



pr< >l>l<-nif. in liadn hi [i!i\ sics( lor a review si e [ijj]). The scniilcpt uiiic 1 Un-ay D A //' is studied usiii« the Q( 'D 



sum rules with three point correlation function 



in [3 



Then, the semileptonic decays D + — * K Q *e + v e 



1G] 
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D ^ tt Iv [17], D p We [18| and B — ► D{D*)h 



191 ] are studied in the same framework. In the present 



work we study the semileptonic decay of B s meson to positive parity D s j(2460) meson, i.e, B s — ► D s j(2460)^, 
within QCD sum rules method. Note that, in 2(J, the decay £? s — > Z> so (2317)ii^ has been studied using the 
QCD sum rules. 

The paper is organized as follows: In section II the sum rules for the transition form factors are calculated; 
section III is devoted to the numerical analysis, discussion and our conclusions. 

II. SUM RULES FOR THE B s -> D aJ (2460)iu TRANSITION FORM FACTORS 

The B s — ► D s j transition proceeds by the 6 -> c transition at the quark level. The matrix element for the 
quark level process can be written as: 

M q = ^V cb v 7/t (l - 75 )Z c 7m(1 ~ 75)^ (1) 

In order to obtain the matrix elements for B s — ► D s j{2^QQ)iv decay, we need to sandwich Eq. fl} between 
initial and final meson states. So, the amplitude of the B s — > D s j(2AQ0)iv decay can be written as: 

M = ^V cb V 7/i (l - <D sJ \ c 7/i (l - 75 )6 I B s > (2) 

The main problem is the calculation of the matrix element < D s j | C7„(l — 75)6 | B s > appearing in Eq. j[2]). 
Both vector and axial vector part of c 7^(1 — 75)6 contribute to the matrix element considered above. From 
Lorentz invariance and parity considerations, this matrix element can be parameterized in terms of the form 
factors in the following way: 

< D s .j(p>, s) I c 7M 75& I B s {p) >= 7 ^ g2) E ^ w P y (3) 

(IUB, + m D sJ ) 

< D sJ (p',e) I c 7m 6 I B s (p) > = i [fo{q 2 )(m Bs + mu.jjej 



(e*p)Pf, + ■ ?(e*p)qn 



(4) 



{m Bs +m DgJ ) (m Ba +m DsJ ) 

where fv(q 2 ), fo(q 2 ), f+{q 2 ) an d f-(q 2 ) are the transition form factors and P M = (p+p')n, q^ = (p — p')^. In 
all following discussions, for customary, we will use following redefinitions: 

fvtf) = , fV j_ q 2 , > foil") = Mq 2 )(m Bs +m DsJ ) 

(m Bs + m DsJ ) 

f + (i 2 ) - , u f ] s , fM 2 ) - . f i q2) , (5) 

(m Bs + m DsJ ) (m Bs + m DsJ ) 

For the calculation of these form factors, QCD sum rules method will be employed. We start by considering 
the following correlator: 

n^ A (pVW) = i 2 J d^ye-^e^'y < | T[J vDsJ (y) J^ A (0) J Bb (x)} | > (6) 
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where J V D s .,(y) = ~sivlbC, Jb s ( x ) = ^75 s i Ja — C7 M 6 and — 07^756 are the interpolating currents of the 
D s j, B s , vector and axial vector currents respectively. 

To calculate the phenomenological part of the correlator given in Eq. ©, two complete sets of intermediate 
states with the same quantum number as the currents Jd sJ and Jb 3 respectively are inserted. As a result of 
this procedure we get the following representation of the above-mentioned correlator: 

ttV,A/ 2 12 2\ 

< I r Dsj I D sJ (p'r >< Dsjjp'Y I Jl A I B s (p) >< B s (p) \J Bs \0> 
(p' 2 -m 2 Dsj )(p 2 -m% s ) 

(7) 

where • • • represent contributions coming from higher states and continuum. The matrix elements in Eq. ([7]) 
are defined in the standard way as: 

<0 J£ sJ D s j(j>')>=f DsJ m DeJ e", <B s (p) \ J Bs >= -i ± s B ° (8) 

m b + m s 

where fo SJ and fs, are the leptonic decay constants of D s j and B s mesons, respectively. Using Eq. ([3|), Eq. 
(gj and Eq. © and performing summation over the polarization of the D s .j meson, Eq. (|7|) can be written as: 



n„„(p 2 ,p' 2 , q 2 ) — ■ — ~t ~/ jo 2° 7 \f 2 2~T x \f'o9iJ.v + f'+P^Pv 

{nib + m s ) (p' z - mj) sj )(p z - m Bs ) 

+ I'-1\-iVv\ + excited states. 
IT 4 (v 2 v' 2 a 2 ) - -is flD '«/^l^ fD.j m D.j f i , 

excited states. (9) 



In accordance with the QCD sum rules philosophy, H^(p 2 ,p' 2 , q 2 ) can also be calculated from QCD side 
with the help of the operator product expansion(OPE) in the deep Euclidean region p 2 <C (mt + m c ) 2 and 
p' 2 <C (m c + m s ) 2 . The theoretical part of the correlator is calculated by means of OPE, and up to operators 
having dimension d = 6, it is determined by the bare- loop and the power corrections from the operators with 
d = 3, < ipip >, d = 4, m s < ipt/j >, d = 5, to 2 , < tpip > and d — 6, < ipiptpip >. In calculating the d = 6 
operator, vacuum saturation approximation is used to set < tptp-tpip > — < ipip > 2 . In calculating the bare-loop 
contribution, we first write the double dispersion representation for the coefficients of corresponding Lorentz 
structures appearing in the correlation function as: 



jjper _ _ 1 f fads'- — P ^ S [ S ] q - + subtraction terms (10) 
(2n) 2 J (s - p 2 )(s' - p' 2 ) 



The spectral densities Pi(s, s', q 2 ) can be calculated from the usual Feynman integral with the help of Cutkosky 
rules, i.e. by replacing the quark propagators with Dirac delta functions: p 'il m 'i ~ * —2n8{p 2 — m 2 ), which 
implies that all quarks are real. After standard calculations for the corresponding spectral densities we obtain: 

p v (s,s',q 2 ) = N c I (s,s',q 2 )[m s + (m s ~m b )B 1 + (m s +m c )B 2 }, 



where 



International Conference on Hadron Physics, 30 Aug - 3 Sep 2007, Canakkale, Turkey 

pais^s^q 2 ) = NJ (s,s\q 2 )[%(m b -m s )Ai-Am b rn c m, s 

+ 4(m s — mj + m c )m 2 — 2[m s + m c )(A + m 2 s ) 
— 2(m s — mf,)(A' + m s ) + 2m s u] 
p + (s,s',q 2 ) = N c I (s,s',q 2 )[4{m b -m s )(A 2 +A 3 ) + 2(m b -3m s )B 1 

-2(m c + m s )B 2 - 2m s ], 
p^(s,s',q 2 ) = N c I (s,s',q 2 )[4(m b -m s )(A 2 -A 3 )-2(m b + m s )B 1 
+2(m c + m s )B2 + 2m s ] 



\{s,s',q 2 ) = s 2 + s' 2 + q 4 - 2sq 2 ~2s'q 2 -2ss', 

A' = (s' -m 2 c +m 2 s ), 

A = (s-m 2 +m 2 ), 

u = s + s' — q , 



A, = 



- ^-[A' 2 s + 2A'm 2 s + m 4 s s + A 2 s' + 2Am 2 s' 

s',q 2 ) 



2A(,s,,s',g2) 



4m 2 ss' — AA'u — Am 2 u — A'm 2 u — m^u + m 2 u 2 ], 

A 2 = — - „ - \2A' 2 ss' + 4A 'miss' + 2mUs' + 6AV 2 
X 2 (s,s',q 2 ) 1 

+ 12Am 2 s' 2 + 6m^s' 2 - 8m 2 s ss' 2 - 6AAVu 

-6Am 2 s'u - 6A'm 2 s'u - 6m 4 s s'u + A l2 u 2 + 2A'm 2 u 2 

+m\u 2 + 2m 2 s s'u 2 ], 

A 3 = 1 , [4AA / aa / + 4Am 2 ss / + 4A , m 2 ss / + Amtss' 

X 2 (s,s',q 2 )' 

—3A' 2 su — 6A'm 2 su — 3m^su — 3A 2 s'u — 6Am 2 s'u 
-imls'u + Am 2 s ss'u + 2AA'u 2 + 2Am 2 s u 2 + 2A'm 2 s u 2 
+2m A s u 2 - m 2 s u 3 } 



(11) 



(12) 



The subscripts V, and ± correspond to the coefficients of the structures proportional to is llva ^'p la V > 9p.v an d 
\{p^Pu ^p'uPv) respectively. In Eq. (fTTj) 7V C = 3 is the number of colors. 
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The integration region for the perturbative contribution in Eq. (|10[) is determined from the condition 
that arguments of the three S functions must vanish simultaneously. The physical region in s and s' plane is 
described by the following inequalities: 



_ 1 < 2ss' + (s + s' - g 2 )(m 2 - s - m 2 ) + (m 2 - m 2 c )2s 

A 1 /2( m 2 jS;TO 2 )A l/2( S)S / )g 2) -+ K ) 

For the contribution of power corrections, i.e. the contributions of operators with dimensions d = 3, 4 and 5, 
we obtain the following results: 

fP + if + if -±,<*>-^<*> { ^ + %a 

m 2 ,2m 2 m 2 + m 2 — q 2 2m?. 

+^ < * s > + r , 2 ; 2 + Tvf] 

m 2 , r 3m 2 3m? 2 

— ^<^>[^ + ^ + ^ 

2m 2 + 2m 2 + m b m c — 2g 2 

' r /2 r 2 J 

/(3) , /(4) /(5) _ (m b - m c ) 2 - q 2 
Jo t Jo -t- Jo - 



m s ^ _ ^ 2m b m 2 + m c m\ + mj; — m c g 2 
m c + m b 2m c m 2 - m| - m b m 2 + m b9 2 



rn 2 f -16m b m;5 + 8mim( + 8to* — 8m 2 o 2 

< ss > J 



+ 



16 [ r' 3 r 

-16m|m c + 8m^ + 8m 2 m 2 — 8m 2 



4m 2 — 8m b m c + 4m 2 — 4g 2 



4m 2 — 8m b m c + 4m? — 4q 



6 

r'r 



H — [ _8m & m c - 8m b m 3 + 8m b m c q 2 + 4m b 



+ 8m 2 m 2 + 4m* - 8m 2 9 2 - 8m 2 g 2 + 4q 4 ] } 
_ml _ r 3m 2 (mg + m 2 -2m b m c -g 2 ) 
12 < SS > [ r'3r 

2 m 2 + m 2 - 2m b m c - q 2 
+om h — = 



-3m b m c (m 2 + ml — q 2 ) + 2(m 2 + m 2 — g 2 ) 2 — 2m 2 m b 



3m c (m c - mi ) + 2(m 2 - g 2 ) 



3m b (-3m c + m b ) + 4(m 2 - q 2 ) 2_ 



/(3) , /(4) . /(5) 1 ^ _ ^ . m s _ -m c m b 

2 _ r 16m 2 16m 2 16 —8m 2 — 8m 2 + 



m 

32 <W> 1- 



-.2 
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ml r 3m 2 3m? 2 

^ 2m| + 2m 2 + m\,m c — 2q 2 ^ 

,'(3) . ,'(4) . /(5) 1 ^_ ^ m s -m c m b 

f- + /- +/- = ^7 < ss > — — < ss > {—^ + —} 

m 2 r 16m? 16m? 16 8m? + 8m? — 8g 2 

1 < ss > r 1 6 1 1 6 1 



m? r 3m 2 3m? 

2m^ + 2m 2 + m b m c — 2g 2 



(14) 



where r = p 2 — m 2 , r' = p' 2 — m 2 . We would like to note that the contributions of operators with d = 6 are 
also calculated. Numerically their contributions to the corresponding sum rules turned out to be very small 
and therefore we did not present their explicit expressions. Note also that, in the present work we neglect the 
a s corrections to the bare loop. For consistency, we also neglect a s corrections in determination of the lcptonic 
decay constants fs a and /d sJ - 

The QCD sum rules for the form factors f v , f' Q , f' + and f'_ is obtained by equating the phenomenological 
expression given in Eq. ([9]) and the OPE expression given by Eqs. (|1HI14[) and applying double Borel transfor- 
mations with respect to the variables p 2 and p' 2 (p 2 — > M 2 ,p' 2 — > M 2 ) in order to suppress the contributions 
of higher states and continuum: 

, ^ = (m b + m s ) 1 e m%JMl + ml 3j /Ml 

fB s m 2 B fD 3j m DsJ 



1 



(27T) 



2^ 



ds / ds' Pl { S , S ',q 2 )e- s ' M "- s 'l M - 

(mb+m s ) 2 J (m c +Ti s ) 2 

+S(if + /i 4) + /i 5) )] 



(15) 



where i = V, and ±, and B denotes the double Borel transformation operator. In Eq. (fl5|) . in order to subtract 
the contributions of the higher states and the continuum, quark-hadron duality assumption is used, i.e. it is 
assumed that 

phigherstates^ ^ = p OPE ^ _ _ g /) (16) 

In calculations the following rule for double Borel transformations is used: 



T(m) r(n) 



m 2 b /M 2 m 2 JM' 2 



(17) 
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III. NUMERICAL ANALYSIS 

In this section we present our numerical analysis for the form factors fv{q 2 ), foil 2 ), f+(q 2 ) an d f-(q 2 )- From 
sum rule expressions of these form factors we see that the condensates, leptonic decay constants of B s and D s j 
mesons, continuum thresholds sq and s' and Borel parameters M 2 and M| are the main input parameters. In 
further numerical analysis we choose the value of the condensates at a fixed renormalization scale of about 1 GeV. 
The values of the condensates are|21| : < Ifip | M=1 GeV>= -(240 ±10 MeV) 3 , < ss >= (0.8±0.2) < ip ip > and 



ml = 0.8 GeV 2 . The quark masses are taken to be m c (p = m c ) = 1.275 ± 0.015 GeV, m s (l GeV) ~ 142 MeV 



22J andm 6 = (4.7±0.1) GeV 



2 11 ] also the mesons masses are taken to be mu sJ = 2.46 GeV and rriB s = 5.3 GeV . 
For the values of the leptonic decay constants of B s and D s j mesons we use the results obtained from two-point 
QCD analysis: f Bs = 209 ± 38 MeV Q and f DsJ = 225 ± 25 MevQ- The threshold parameters sn and s' 



12j. The 



are also determined from the two-point QCD sum rules: sq = (35 ± 2) GeV 2 [13[ and s' = 9 GeV 2 
Borel parameters M 2 and Mf are auxiliary quantities and therefore the results of physical quantities should not 
depend on them. In QCD sum rule method, OPE is truncated at some finite order, leaving a residual dependence 
on the Borel parameters. For this reason, working regions for the Borel parameters should be chosen such that 
in these regions form factors are practically independent of them. The working regions for the Borel parameters 
Ml and M| can be determined by requiring that, on the one side, the continuum contribution should be small, 
and on the other side, the contribution of the operator with the highest dimension should be small. As a result 
of the above-mentioned requirements, the working regions are determined to be 10 GeV 2 < M\ < 20 GeV 2 and 
4 GeV 2 < M 2 < 10 GeV 2 . 

In order to estimate the width of B s — > D s jlv it is necessary to know the q 2 dependence of the form factors 
fv(q 2 ), fo(q 2 ), f+(q 2 ) an d f-(q 2 ) in the whole physical region mf < q 2 < (raj - wd sJ ) 2 - The q 2 dependence 
of the form factors can be calculated from QCD sum rules (for details, see [3, For extracting the q 2 

dependence of the form factors from QCD sum rules we should consider a range q 2 where the correlator function 
can reliably be calculated. For this purpose we have to stay approximately 1 GeV 2 below the perturbative cut, 
i.e., up to q 2 = 8 GeV 2 . In order to extend our results to the full physical region, we look for parameterization 
of the form factors in such a way that in the region < q 2 < 8 GeV 2 , this parameterization coincides with 
the sum rules prediction. The dependence of form factors fv(q 2 ), fo(q 2 ), f+(q 2 ) and f-(q 2 ) on q 2 are given in 
FigsfTJ [21 [3] and 0] respectively. Our numerical calculations shows that the best parameterization of the form 
factors with respect to q 2 are as follows: 



Mq 2 ) = ^- (18) 

1 + aq + f3q 2 + 7<f + A<? 4 

where q — q 2 /m 2 B . The values of the parameters fi(0), a, (3, 7, and A are given in the Table 1. 
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f(o) 


OL 





1 


A 


fv 


1.18 


-1.87 


-1.88 


-2.41 


3.34 


fo 


0.076 


1.85 


0.89 


19.0 


-79.3 


/+ 


0.13 


-7.14 


11.6 


21.3 


-59.8 


/- 


-0.26 


-4.11 


-3.27 


15.2 


18.6 



TABLE I: Parameters appearing in the form factors of the Bs — ► D a j (2460)£v 
decay in a four-parameter fit, for M\ = 15 GeV 2 , M| = 6 GeV 2 

For B s — > D s j(2460)Zi/ decay it is also possible to determine the polarization of the Z? s j (2460) meson. For 
this aim we determine the asymmetry parameter a, characterizing the polarization of the D s j(2460) meson, as 



^T L /dq 2 
a = 2— — r-rrr — 1 



(19) 



dY T /dq 2 

where dT L / dq 2 and dT T / dq 2 are differential widths of the decay to the states with longitudinal and transversal 
polarized D s j (2460) meson. After some calculations for differential decay rates dTL/dq 2 and dTr/dq 2 we get 



dT T 



7\G%\ V cb | 2 {(2A + Bq 2 )[\ f v | 2 (4m| s | 7 ?)+ I & | 2 ]} 



dq 2 87r 4 m B 2 

16 J mB2 ^\ G F\ V ^ 2 {(^ + Bq 2 ) [\ fv I 2 (4m| s | ^ 
+m| 



2 I P | 2 / 2 2 2\\ i I rl |2 

-(m B -ran, -<?))+ I /o I 



- 1 /; i 2 ^%^(2mi s + 2^ sj - {)- 1 /i i 2 

l/tn l/tn 



™ 2 Bs 1 7 12 



-2 ' - ' (fle(/o/+ + + (m Bs - m^J/;/!*)) 



2 I 7? |2 



B s I P I h J-/ |2 , / 2 2 \2 I rl |2 , 4 I rl |2 

~ 2i? 7^2 LI JO I + (™B S - m D sJ ) I /+ I +<? I /- I 

-D sJ 

+2( m|, - m 2 DaJ )Re{f'J+) + 2g 2 /^:* + 2? 2 (m| s - ro^ife^/l*)] } 



(20) 



(21) 



where 



7 



A 
B 

h 



AV2( T 



,9 2 ) 



2m Bs 



— (g 2 -m 2 )( g 2 + 2m 2 )/ 



12 g 2 
1 
6~q 



= 77(1 --f) 



(22) 
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The dependence of the asymmetry parameter a on q 2 is shown in Fig. O From this figure we see that asymmetry 
parameter a varies between -0.3 and 0.3 when q 2 lies in the region m 2 < q 2 < 6 GeV 2 . An interesting observation 
is that around q 2 = 5.2 GeV 2 the asymmetry parameter changes sign. Therefore measurement of the polarization 
asymmetry parameter a at fixed values of q 2 and determination of its sign can give unambiguous information 
about quark structure of D s j meson. 

At the end of this section we would like to present the value of the branching ratio of this decay. Taking 
into account the q 2 dependence of the form factors and performing integration over q 2 in the limit mf < q 2 < 
(niB s — mD SJ ) 2 and using the total life-time tb s = 1-46 x 10~ 12 s [27[ we get for the branching ratio 

B(B S -» D sJ (2Am)tv) ~ 4.9 x 10~ 3 (23) 

which can be easily measurable at LHC. 

In conclusion, the semileptonic B s — > D s j(2AQQ)£v decay is investigated in QCD sum rule method. The q 2 
dependence of the transition form factors are evaluated. The dependence of the asymmetry parameter a on q 2 
is investigated and the branching ratio is estimated to be measurably large at LHC. 
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FIG. 1: The depend* 
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FIG. 2: The dependence of f on q 2 at M 2 = 15 GeV 2 , Aff = 6 GeV 2 , s = 35 GeV 2 and s' = 9 Gel/ 2 . 




FIG. 5: The dependence of a on q 2 . 



